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Abstract 

The solution = ('!9f)^>Q of a class of linear stochastic partial differen- 
tial equations is approximated using Clark's robust representation approach 
([1], [2]). The ensuing approximations are shown to coincide with the time 
marginals of solutions of a certain McKean-Vlasov type equation. We prove 
existence and uniqueness of the solution of the McKean-Vlasov equation. The 
result leads to a representation of as a limit of empirical distributions of 
systems of equally weighted particles. In particular, the solution of the Za- 
kai equation and that of the Kushner-Stratonovitch equation (the two main 
equations of nonlinear filtering) are shown to be approximated the empirical 
distribution of systems of particles that have equal weights (unlike those pre- 
sented in [15] and [16]) and do not require additional correction procedures 
(such as those introduced in [5], [6], [10], etc). 
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1 Introduction 



Let (r2,jF, P) be a probability space on which we have defined an m— dimensional 
Wiener process W = (VF*)^!- Let M{W^) be the space of finite Borel measures 
defined on the (i— dimensional Euclidean space M*^ and {'&t)t>o be an M(]R'^)-valued 
stochastic process satisfying the following linear stochastic partial differential equa- 
tion 

ft m ft 

M^) = M^) + / {as^ + Ls^) ds + y^ iPi^) dWi, (1) 

^0 Jo 

where L — (LJ^^q , Lg : C6(M'^) Cb{M.^) is the second order elliptic differential 
operator 

d d 

^.V^ = 2 E + E ^sd.,"^, (2) 

i,j=l 1=1 

and 9? is a function in the domain of L, (/? G f^^^Q P (Lg). For simplicity, we will 
assume that t?o is a probability measure. If a = 0, then (1) is called the Zakai 
or the Duncan-Mortensen-Zakai equation (cf [11], [21], [27]). The Zakai equation 
has been studied extensively over the last 40 years because of its importance in 
non- linear filtering theory (see [20], [22]): In non- linear filtering, -dt, the normalized 
form of 'dt, is the conditional distribution of a (non-homogeneous) Markov process 
with infinitesimal generator L given the observation process W which satisfies the 
evolution equation 

Wt^ f /3s{^t)ds + Vt, (3) 

In (3), V is an m— dimensional Wiener process independent of ^. Within the filtering 

problem W is not a Brownian motion as it is assumed in the current set-up. How- 
ever, it becomes a Brownian motion after a suitable change of measure (Girsanov 
transformation). One can prove that satisfies the following (non- linear) stochastic 
partial differential equation 

Jo 

+ E / (f^sv) - (Pi) iv) {dwi - (Pi) ds) , (4) 

J=i 
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If q; = 0, then equation (4) is called the Fujisaki-Kallianpur-Kunita or Kushner- 
Stratonovitch equation (cf. [12], [17]). In general, neither (1) nor (4) have explicit 
solutions, though one can approximate them by numerical means. As expected, 
there is a wide variety of methods to do this (see, for example, [3] and the references 
therein). Among them, particle methods seem to be in many cases quite effective. 
The starting point for a particle approximation of -dt is the following Feynman-Kac 
representation 

^t{v) = E[ip{Xt)A{X)\W], (5) 

where X = {Xt)^^^ is a Markov process with infinitesimal generator L independent 
of W with initial distribution and {X) is defined as 

A {X) = exp a, {X,) ds + J2 (/^ i^s) dW^ - ^ Pi {X,f ds) j , t > 0. 

Since in (5) X is independent of W, one can use the Monte Carlo method to compute 
That is, i&t has the representation 

1 " 

i9t= lim -y2A{X^)d^., (6) 

1=1 

where X\ i > are independent realizations of X. In other words, {XI, At {X'^)}°^^ 
satisfies the following system of SDEs 

r dXl =ht{Xl)dt + at{Xl)dBl 

\ dA {X^) = A {X^) at {XI) dt + Er=i At {X,) (51 {XI) dW} ^' > 

where S*, i > 0, are mutually independent rf— dimensional Brownian motions inde- 
pendent of W and = {o's)'ij=i is chosen to satisfy = crs(^J (o'J is the transpose 
of as). As demonstrated in [15] and [16], representations of the type (5) and (6) 
hold true for a far wider class of SPDEs than the one described by (1). 

However, the convergence in (6) is very slow. That is because the variance of 
the weights At (X*) , z > 0, increases exponentially fast with time. The effect is 
that most of weights decrease to with only a few becoming very large. In order 
to offset this outcome, a wealth of methods have been proposed. In filtering theory, 
the generic name for such a method is that of a particle filter ([4], [5], [6], [7], 
[8], [10], etc.). The standard remedy is to introduce an additional procedure that 
removes particles with small weights and adds additional particles in places where 
the existing one have large weights. Put differently, one applies at certain times a 
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branching procedure by which, each particle will be replaced by a random number of 
"offsprings" with a mean proportional with its corresponding weight. This branching 
procedure is a double edge sword: applying it too often may actually decrease the 
rate of convergence (cf. [9]). 

In this paper, we seek a different remedy to the slow convergence of the Monte 
Carlo method. Heuristically we seek to keep the weights of the particles equal 
without introducing an additional procedure but only by amending the motion of 
the particles in a way that will take into account the state of the entire system. This 
will be achieved in two steps: 

First we will show that "dt and its normahsed version -dt admit a robust version 
following the approach first introduced by Clark ([1], see also [2]). By this we mean 
that dt and, respectively, dt depends continuously on the generating Brownian path 
s — * Ws {uj) , uj E Vt. Hence if we consider an approximating sequence of paths 
s — >• Wg iuj) that will converge to s — >• Wg iuj) as 5 tends to 0, then the corresponding 
measure valued processes "§1 and '&1 will converge to 'dt and, respectively, "dt with a 
certain rate of convergence r(t, W. [to) , 5). 

The second step is to prove that ■dl has the asymptotic representation 



= hm where t^^'' 



1 " 



and X*, i > Oare non-linear diffusions satisfying the non- linear SDE 

dXi'' = hi (^f , Xi'') dt + at [Xl'') dBl (9) 

In (9), "dl is the conditional distribution of X]'^ , are mutually independent 
Brownian motions and the function will depend intrinsically on the chosen path 
s 1^5(0;). We will show that '^f'^ — df"' and that, by suitably choosing the 
parameter 5 — 5{W. {u>) ,n), we will have 

1 " 

■dt — lim where -df^ — — } S i,s(w.(aj),n). 

n— >oo n ' -^t 

i=l 

Hence t^t has an asymptotic representation involving particles with equal weights. 
Finally in order to obtain the corresponding asymptotic representation for i.9twe 
need to "unnormalize" d", i.e., to attach to each particle a weight a" (the same one) 
which converges to as n tends to infinity. For example, one could choose 

< = exp n * ds+Y, i^f^ ^'r {P's) d {wty m'ds 
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Hence, we obtain the asymptotic representation i?t which uses systems of particles 
with equal weights. 

"lit — lim ~ / ^ ^^,&{w.^w),■n■). 

n— >oo n ^ — ' t 

i=l 

We remark that we are not aware of a similar result even for the simpler case 
when (1) has no stochastic term, i.e., (1) is a second order elliptic PDE with a zero- 
order term. In this case the first step described above is not required: One obtains 
directly an asymptotic representation of the solution 

1 " 

■dt = lim - V 5 V. 

i=l 

where the X*, i > are non-hnear diffusions satisfying a non-hnear SDE of the type 
(9). 



2 The robust representation of '&i 



In this section we introduce the robust representation formula for 'dt and, respec- 
tively, "Qt- By robustness here we mean that the dependence of the generating 
Brownian path t ^ Wt (cu) is continuous. The formula is similar to, and inspired 
by, the robust version of the integral representation formula of nonlinear filtering as 
presented in [1] and [2]. At a formal level, the formula is derived by a process of 
integration-by-parts applied to the stochastic integrals that appear in the Feynman- 
Kac representation (5). The rigorous justification of the formula is identical with 
that of the corresponding result in [2] and for this reason we omit it here. 

For the existence of the robust representation of i^t and, respectively, of -^t we 
follow Theorems 1 and 2 in [2]. For this we need to impose the following conditions: 



X 



^) , hence 



RRl: For all j = 1, c? we assume that (s, x) — >■ (3^ 
s Pi (Xg) is a semimartingale (we denote by (X)™" its martingale part and by 
l3^ (X)"^^ its finite variation part) such that, for all A; > and alH > 



E 



< oo, E 



< oo. 



RR2: For alH > we have 

E 



exp 2 I as (Xg) ds 





< oo. 
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Let y. — {us, s > 0} e C^a ([0,oo)) be a continuous path R'^ and define = 
(Of )^>Qto be the following M (R'^)-valued stochastic process {if is a bounded Borel 
measurable function) 



E 



as{X,)ds + J2pi{X,)yi 

- E (?/.) + 2 I i^liXsfds^ 



where X = {Xt)^^^ is a Markov process with infinitesimal generator L indepen- 
dent of W with initial distribution ??o and P (y.) is a version of the stochastic inte- 
gral yid/Si (Xs). Let B^' = (Gf j^^pbe normalized version of the M (M-^) -valued 
stochastic = (Bf )^ 



't>0 



* ©r(i)' 

where 1 is the constant function l[x) = 1 for all a; G M'^ (0^' (1) is the mass process 
associated to ). Following Theorem 1 in [2], under the conditions RRl and RR2, 
and©^' depend continuously on the path y.. Moreover, for any R > 0, there 
exist two constants K — K [R, t) and K — K [R, t) such that 



(10) 



sup 



for any bounded measurable function Lp 
any two paths y.^ such that 1 lt , 1 lt < R (where 
maxi<j<dmaXsg[o,t] |«*|). Furthermore, if we use the norm 
signed measures 



(x)| < oo) and for 
1^ is defined as \\a\\^ = 
■IL„ on the set of finite 



l/^IL= sup |/u((/j)|, 

{<^eC6(Md)|||<^||<i} 



then, from (10), we deduce that 

<K\W-y''\ 



Of - ef 



ef - ef 



<K\\y'-y' 



(11) 



Following Theorem 2 in [2], 0*^ and0^ are the robust versions of i} and, re- 
spectively. More precisely for all t > 0, 







w. 



and e^ = '&t-i P-almost surely 



(12) 
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and the null set can be taken to be independent of t since all processes involved are 
time-continuous. 

The robust representation result also enables us to fix the generating Brownian 
path s — > Ws (io) and proceed to approximate the corresponding and, implicitly, 
for this fixed, but arbitrary, path s — > Wg {u). 

We will replace the generating Brownian path i — > (a;) with a smoother one. 
For example, we choose a partition of the time interval [0, oo) of the form 

= to <ti < .... <ti < ... 

where ti — iS, i >0 and define the following piecewise linear path t {uj) 

W! {u) = W,, (^) + (^) - (^) _ , for t e [t,, t,^^). 

Let ^\ = ef ' and ^ = ef Then, from (11), it follows that. 



lim ^\ = and lim 



5^0 



5^0 



where the above convergence is in the weak topology on M (IR'^) . Moreover, we have 
the following rates of convergence: there exist two constants K = K {wt (cu) , t) and 
K — K [vDt {oj) , t) where Wt {oj) — supggjQ^^.,.^] maxi<j<d \ Wl {uj)\ such that 

1 1^,^ - t?J I „ < Kwl (a;) , 1 1^,^ - M < Kwl {u) 



and Wf. {uj) 



maxj=o,[t5] ^^'9te[u,u+,_\ maxi<i<d 



Wl{uj)-Wl {u) 



Moreover since 



s — >• Ws (a;) is continuous, we can choose 5 = 6 {uj,n) , so that {u) < -. Hence 



In (13), the norm 



1?; 



\M. 



5{u!,n) 



< 



M 



< 



M 

^ is defined as 



— V" 
K 



< 



(13) 



k=0 



where ipo — 1 and {ipk)k>o are the elements of a dense set A4 E Cb 
will be used to prove convergence in the next step of the construction. 

^Both norms \ \-\\^ and .|H|^ induce the weak topology on M 



(14) 
This norm 
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It is easy to see that -^f satisfies the following hnear PDE (written in weak form) 
i^tiif) = Mv) + f < + Ls^) ds (15) 

where 

whilst "dl satisfies the following nonlinear PDE (again, written in weak form) 

- ^o(^) + f («f ^ + Ls^) - W) {V) ds (17) 

Note that, from (15), ■dl{l) = exp J*^l (af) ds^ (since ^^^{1) = 1). 

We introduce now additional conditions on the coefficients of the SPDE (1) to 
insure it has a unique sohition and that the solution has a density with respect to 
the Lebesgue measure and its density is strictly positive everywhere. 

EU: We will assume that the functions ag, Pi, j — l,...,m, are bounded (with 
a bound independent of s). al?,i,j = l,...,d, h\, i = l,...,d are bounded (with a 
bound independent of s) and Lipschitz (with a Lipschitz constant independent of s). 
We also assume that the operator L is uniformly elliptic. That is there is a constant 
b such that for all times s > and vectors ^ = (6)iLi ^ I^*^) we have 

where (■, ■) is the standard inner product on M.^. Further, we will assume that t^q has 
finite second moment. 

PD: We will assume that -do is absolutely continuous with respect to the Lebesgue 
measure and its density is strictly positive everywhere. 

Assuming condition EU ensures that the SPDE (1) (see for example Jie&Kurtz) 
and the PDEs (15) and (17) have unique solutions. Further, the system of equations 
(7) and, in particular, the SDE satisfied by the processes (which appear in 
the asymptotic representation (6) of i^t ) has a unique solution. In particular this 
ensures the existence of a Markov process X = {Xt)^^^ with infinitesimal generator 
L and initial distribution 'd^. Assuming EU+PD guarantees that, for any t > the 
distribution of Xt is absolutely continuous with respect to the Lebesgue measure 
and its density is strictly positive everywhere. 
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Let us note that, from the Feynman-Kac representation (5), -dt (and hence its 
normaUzed version "i^^) has the same support as the distribution of X^. Similarly, 



'dl has the following representation 



(^) = E 



V? {Xt) exp 



[Xs) ds 



(18) 



So, both ■&1 and have the same support as the distribution of Xt^ too. It follows 
from the above that 1?^, ■^j, and arc all absolutely continuous with respect to the 
Lebesgue measure and their density is strictly positive everywhere. In the following 
we will denote by a; — >• (x) and x ^ 'dl {x) the density of 'dt and respectively 
■§1 with respect to the Lebesgue measure. 



3 The non-linear process 

We introduce next the non- linear SDEs satisfied by the non-linear diffusions X*'*^, 
i > appearing in the asymptotic representation (8) of 'd^. For this we need to 
define first the drift coefficient of the non-linear SDE (9). In the following, we fix 
the generating Brownian path t ^ Wf {u) . The coefficient af obviously depends on 
the path, however we will not make the dependence explicit. For this fixed path, 
the non-linear diffusions X^'^ will be defined on a new probability space (^Q, T ^ . 

For arbitrary f eB (M'^) and e P (M"') , define f^eB (M'^) to be the function 
fi^^f-lj,(f) and let A J = {AifY-^^ be the vector function 



If {y - x) jy) 

JR-I \\x — y\\ 



y\ 

where Ud is the surface area of the d— dimensional sphere S^-i- li d — 1, then A^/ 
has several equivalent representions representation 

A^/ {x) = (/l[x,oo)) - ^ (/) ^ (l^oo)) = ^ (/''l[x,oo)) 

= 1^ (l(-oo,x)) /"(/)-/" (/l(-oo,x)) = -/W {rM-oo,x)) (19) 



When d — 1 the function A^/ is well defined for arbitrary f E B (R") and /x e 
V (M'^). This is not the case when d > 2. Let fi E V {R'^) , d > 2 be a probability 
measure absolutely continuous with respect to the Lebesgue measure. We say that 
/X is a p—good probability measure if its density x ^ d^ {x) with respect to the 
Lebesgue measure is locally bounded and p-integrable, where 1 < p < d. The 
following proposition gives a necessary condition which insures that A^/ is well 
defined: 
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Proposition 1 If fj, E V (M"j ,d>2 is a p-good probability measure, then K^f is 
well defined. 



Proof. We need to prove that the function 

y-x 



y 



\x - y\ 



ifiy)-f^{f)) 



is intcgrable with respect to fj, for all a; > 0. Let ||/|| be the L°°-norm of / and 



f^\B{x,l) 



be the L°°-norm of the function x ^ d^ {x) restricted to B{x,l) , the 



Dall of radius 1 and centre x. Then 



1 



\x - y\ 



i^ (dy) < 



! I 

/ II 

Jb{x,i) \\x-y\\ 



dy 



+ 



B{x,ir \\x-y\ 



ZTt^idy), (20) 



where -B(a;, 1)^ = M.^\B (x, 1) . By using polar coordinates it is easy to check that 
the first integral is equal to Ud and hence it is finite. Using Holder's inequality and 
polar coordinates, we get that 



B{x,ir IF - 2/1 



7[i^(dy) < 



B{x,iy' \ \x — y 



dy 



B{x,ir 



{dMYdy 



< (iOd r-M^'^-^^drV 11^ 



and since (g — 1) {d — 1) — ^ > 1 we get that r~(^~^)^'^~^)<ir is finite. Hence 
the second integral on the right hand side of (20) is finite and so A^/ is well defined 
and is a bounded function. ■ 



We now define the first order differential operators 



d^{x) 



In the following, we will denote by 1^ to be the Lebesgue measure on 
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Proposition 2 

i. If ^pe Cl (R), then 



a. If d > 2, E cl (R'^) and fj, E V [W') ,d>2 is a p-good probability measure 
(1 < p < d), then 

Proof. 

i. Using the representation (19) we have that hm|3,|^oo i^) — 0, hence by inte- 
gration by parts 

/ if {x) f^ {x) n {dx) = / ip'Au_f {x) dx= Lj^ip {x) dx = U^^p {x) n {dx) . 
Jr Jr Jr. Jr 

ii. For (p E Cl (M'^) the following identity holds true (see, for example, [26] page 12) 



Jr 

We will show that the function 

{x, y) eR"^ xW^ — > 



1 f {y - x,Vip{x)) 



{y-x,V(p(x)) 



dx. 



\\x-y\\'^ 

is integrable with respect to l2d- Hence by Fubini's theorem 

f^ivf') = / ^{y)r(y)d,(y)dy 

1 f {y - x,Vip{x)) 



r {y) d, {y) 



i_ f {y~x)f^{y) 

\ 



^ dxf (y) d^ (y) dy 
d (y) dy, (x) 



\\x-y\\ 

(A^/(x),V^(x))= / Ll^{x)dx. 

jR'i 



First let us observe that 

{y-x,Vip (x)) 



\x-y\\' 



F (y) d, (y) 



< 2 



\\^-y\\ 
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Let now h be the following Riesz type operation 



Jmfi \ \x — v\\ 



x-y\ 

for ip e C'' (R'^). Theorem 1.2.1 from [26], page 12, states that hip is g-integrable 
for any q such that 

111, , , 

- = : where 1 < p < a. 

q p a 

So Iiip is g-intcgrable for any q such that < - < 1 — |. Hence Iiip is g-integrable 
for q such that ^ = 1 — K So, by Fubini's theorem (for non- negative functions) and 
Holder's inequality, 

JR'ixR'i \\x — y\\ jRd \Jud \\x — y\\ J 

hi^ (y) d„ (y) dy 



< (I hi;{yYdy\(l d,{yYdy\ 



where [x) = \\^(p {x)\\. Hence our claim. ■ 

We are now ready to define the coefficients fe^of the equation (9). Let /x be a 
probability measure that satisfies the following conditions: 

1. is absolutely continuous with respect to the Lebesgue measure and its 
density x (x) with respect to the Lebesgue measure is a strictly positive 
function in Cl{M!^). 

2. If (i > 2, then there exists p such that is a p—good probability measure. 
Then the following coefficients are well defined 

(5, x) hi x) ^h{x)-r ^^^^J = 1, m. (21) 

Let X = {Xt)^^Q be a continuous (i— dimensional stochastic process. In the 
following we will denote by 1?^ the distribution of X^. We say that X is a good 
process if the following three conditions are satisfied: 

• is absolutely continuous with respect to the Lebesgue measure and its den- 
sity x — > (x) is positive for all x e W^. 
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• The function (t, x) e [0, oo] x R*^ — >■ i^t (x) is continuous. 

• If (i > 2, then there exists p such that i9t is a p—good probabihty measure for 
all t > 0. 

We have the following proposition 

Proposition 3 Let (q,^,P^ be a probability space on which there exists a good 
process X = which satisfies equation (9). That is 

Xf = X^ + hi (^f , X'^) ds + a, {X',) dB, 



where the coefficients bg are those specified in (21), B = {Bt)^^^ be a d- dimensional 
Brownian motion and Xq has distribution do. Then i?^ will be equal to d^. 

Proof. Let us apply Ito's formula to the equation (9) for (p e (IR'^). We get that 




which yields, by taking expectation and applying Proposition 2, that i?* satisfies 

= ¥M + /* di {aiv + L,<p) - di {ai) ^i i^) ds 
Jo 

So satisfied the PDE (17) and, using the uniqueness of the solution of (17), it 
follows that d^ ^'d^.m 

4 Uniqueness of the Solution 

In the following we will prove the uniqueness of a solution of (9) in the class of good 
of processes as defined in the previous section. 

Theorem 4 There exists at most one solution of (9) which is a good process. 
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Proof. First, we note that if — [Xf)^^^ is a solution of (9) then the distribution 

of Xf satisfies (17) and therefore is uniquely determined and equal to -dl- Therefore 
we only need to justify the uniqueness of the solution of 

dr/ = hi (dl Xi:") dt + at (X;) dBl (22) 

which is obtained from (9) by replacing with §1. Hence (22) is an ordinary SDE 
whose coefficients are locally Lipschitz. In particular x is locally Lipschitz. 

Hence, for example by Theorem 3.1 page 164 in [13], there exists a stopping time C, 
such that equation (9) has a unique solution in the interval [0, Q and, on the event, 
C < cxD we have 

limsup I \xf \ I = oo . 

We want to prove next that the event C, < oo has null probability (hence the solution 
is unique for all t > 0). To do this, it is enough to show that, for alH > we have 
I I < oo P— almost surely. This fact implies, by the continuity of the trajectories, 
that sup[o^t] I \ \ I < oo almost surely and hence that C ^ ^ almost surely. For this 
it suffices to prove that 



lim P(||xf|| > A;) = 0. (23) 



We have 



where 



I 1 if ||x|| > A; 

(x) ^ Inm therwise^ ' ' \ ( Nhf ) ^ < ' '^1 ' < ^ 

' 11x11 <'= 



2 



So (23) is implied by 



lim<(^fe) = 0. (24) 

k—^oo 



Now since the coefficients of the operator Lg are all uniformly bounded it follows 
that (fk £ r\s>o ^ i-^s) a-nd that limfe_^.oo sup^gjQ | |Ls(^fe| | = 0. For arbitrary T > 0, 
define now the function * : [0, T] ^ M+ 

*(t) = sup i}l{^k)- 
se[o,t] 
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Prom (17) one deduces that, for arbitrary T > and t e [0,T] we have 

fT nt 



se[o,T] Jo Jo 



ds, 



where 

^T,a,/3,a; = SUp ||q!,|| 

se[o,T] 



+ -J](sup + sup \\(3i\\ max {Wl^^ {uj) - Wl {co))) 



Hence, by Gronwall's inequahty, 

fT 

ll-^sV^fell / 

se[o,t] \ se[o,Ti 

which imphes that 



sup ^?,^(^fc)<e^W^(4(^^)+ sup / ^\{\)ds 



hm sup ^\{i^k) = 0. (25) 

as hmfc^oo ^o(</'A:) = 0. But (25) imphes (24), which in turn imphes (23) and that 
completes the proof of the Theorem. ■ 



5 Existence of the Solution 

Wc are now ready to complete the last step of the programme. For this we need to 
add one final condition: _^ 

ES: Let {s,x) — > bl'^" (x) be the function defined in (21) where the measure 
II is chosen to be T^f . In other words. 



(26) 



We will assume that (s, x) 



{x) is globally Lipschitz. More precisely we will 



assume that, for any T > 0. there exists a constant Kt such that 



(x) 



(y) 



< Kt \ \x — y\ 



(27) 



for all x.yeR^ and s e [0, T] . 
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Theorem 5 Under the conditions EU+PD+ES, equation (9) has a solution which 
is a good process. 

Proof. We prove the existence of the solution on an arbitrary time interval [0, T] . From 
(27) and the fact that the function 

.e[0,T]^6.(0) + ^^ 

is continuous, hence bounded, it follows (for example, by using Theorem 2.9, p. 289, 
in [14]) that the equation 

dXi'' = b'/' (x;'^) dt + a, [Xl) dBl, (28) 

has a (unique) solution whose distribution satisfies the hnear PDE 

^tiip) = M'P) + £ + L';iip) ds (29) 

which has a unique solution (see ... in [18]). From (15) it follows that -i^f is a solution 
of PDE (29), hence and hence X'^^ is in fact a solution of the nonlinear SDE 
(9).B 

Remark 1 If d = 1 and for all s > 0, the density of-d^ with respect to the Lebesgue 
measure has the form 

(x) = e-^'(^), xeR. 
where is a differentiable convex function, condition ES is satisfied. 

Proof. We have three cases: 

1. If ^{x) = 0, then ^ {^^^^ = {(4)^ which is bounded since af is 
bounded. 

2. If ^ (x) > 0, then using the first part of the representation (19) we get 

i (^) = W)' + [ (^) '4 

Since is convex, iox y > x we have 
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So 

?s 



Jx "'•^ '^•^ J X 



< \\(aiY\ 



since ^ {x) e(--^)i^(-)rfy = 1. 

3. If ^ (x) < 0,then then using the second part of the representation (19) we get 



d fKi)sJ{x) 



dx V e-^'W 



{alY {x) - [alY (y) ^ (x) e^'^^^-^'^^Uy 



and we follow the same steps as in the previous case. ■ 

We complete the section by noting that the above construction works with min- 
imal changes when (the initial Markov process) ^s a reflecting boundary diffusion. 
In this case, the analysis simplifies considerably if the domain is compact. For ex- 
ample, the cumbersome condition ES is replaced by the assumption that ^ has a 
density which is bounded away from 0. We will detail the analysis of this case 
in a forthcoming paper, together with the description of the associated numerical 
algorithm. However, for completeness, we briefly present the results here: 

Following the notation and results in [19], assume that ^ is the solution of a 
stochastic differential equation with reflection along the normal. In other words, ^ is 
the unique solution of the equation 

it = io+ I hs{is)ds+ f <Js{is)dB,-ku (30) 
JO Jo 



where ^ B (0, M)^ for alH > and A; is a bounded variation process 

For general domains D, k is defined as 

l^lt= / Mi.^dD}d\k\s, h= j n{is)d\kl 
Jo Jo 

where n (x) is the unit outward normal to dD at x. 

Obviously the generator L — (ivs)s>o associated to $, has the form (2) for any 
(f in C|- {B (0, M)) , the set of twice differentiable functions defined on B (0, M) 

"^B (0, M) is the ball of center and radius M. 
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with compact sui)port. The operator A^/ is well defined for measures /i with 

support in B (0, M) absolute continuous with respect to the Lebesgue measure and 
with bounded density. The definition of a good process is now slightly simpler: A 
continuous process X = {Xt)^^^ which takes values in B (0, M) is a good process if 

the following two conditions are satisfied ( as before, t?* is the distribution of Xt): 

♦ -dt is absolutely continuous with respect to the Lebesgue measure and its density 

X ^ "dt (x) is positive on B (0, M). 

♦ The function (t, x) G [0, oo]x B (0, M) — > {x) is continuous and bounded. 
Proposition (3) now becomes: 

Proposition 6 Let ^f2,jF, be a probability space on which there exists a good 
process X^ = (^X^)^^^ which satisfies equation 

Xt = X', + £ bl (^f , ) ds + as (Xf ) dB, - kl (31) 
where Xf e S(0,M) for all t > and is a bounded variation process 

pt r-t -^5 

In (31), the coefficients bg are those specified in (21), B = {Bt)^^Q be a d- dimensional 
Brownian motion and Xq has distribution ^q. Then the distribution of X^ will be 
equal to df. 

The uniqueness of the pair (X^, /c^) follows by slight variation of the argument 
for the proof of theorem (4). The existence of the pair [X^, k^^ requires the following 
simpler condition: 

ES': Assume that, ^ = (6)t>o > solution of the SDE (30) is a good process 
and that its density i^t (x) is bounded away from 0. More precisely, for any T > we 
have 

inf (x) > 0. 

t£[0,T]xeB{0,M) 

Then we have the following result: 

Theorem 7 Under the conditions EU+PD+ES', equation (31) has a unique so- 
lution which is a good process. 
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The results follows by proving that ES' imphes ES. Since 1?^ has the represen- 
tation 



ip (e*) exp ^ af (e.) - ^ (af ) 



(32) 



•§1 is absolutely continuous with respect to the Lebesgue measure, too. Further, 
since cuf as defined in (16) is a bounded function on B (0, M), we deduce that there 
exists a positive constant e such that 

{^) > eE iCt)] = e / <fi{x)^t (x) dx 
Jb{o,m) 

Hence the density of is uniformly bounded away from 0. One can also prove that 
the density of df has uniformly bounded first order derivatives. This implies that 



the function (s, x) 



is globally Lipschitz, hence ES is satisfied. 



Remark 2 Numerical methods for equations of type (30) have extensively developed 
(see for example [23], [24] and the references therein). We will adapt these schemes 
in order to produce a numerical method for solving (31). 



6 Rates of convergence and final remarks 

It is obvious that, for any ip E Cf, (M*^) , we have 



E 



< 



m 



so 



E 



< 



where the norm ||-||^ is defined in (14). Prom (13), we have that 



5{L0,n) 



K 

< —= and 
M. \/n 



< ^ 



hence 



E 



5{u),n),n 



M 



< 
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so "d^ ' converges to i^t P- almost surely. We also have almost sure convergence 
if we view and as processes on the product space (f^,^, P) 

(0,^, P) = {nx fi,J^(g)J^,P(g)P) 

on which we 'lift' the processes and from the component spaces. Finally, 

in introduction we chose the asymptotic representation for t?^ to be 

a" ^ 

'&t = lim — > 5„i,5{W.(u),n), (33) 

n-^oo n ^ — ^ t 
i=l 

where 

Jo j^-^ Jo 

-^E / ^t^'^'^'^'^PiT dsj (34) 

Then each of the terms in the formula (34), i.e., J^* (q/^) ds, 

/J (/3^)d(wf"'"^)' and Jo (/3ffds converge P-almost surely to 

/o (ois) ds, Jq iPi) dWi and, respectively, ■dt iPi)"^ ds, hence converges al- 
most surely to (1), so (34) holds true. 
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